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Abstract. We present Giant Metrewave Radio Tele- 
scope (GMRT) HI 21cm emission images of the 
z = 0.009 damped Lyman-a (DLA) absorber to- 
wards the QSO HS 1543-K5921. The DLA has 
been earlier identified as a low surface bright- 
ness galaxy SBS 1543-^593 ( [Reimers fc Hagen 1998 , 
Bowen et al. 2001a) with small (~ 400 pc) impact param- 
eter to the QSO line of sight. The extremely low redshift 
of the absorber allows us to make, for the first time, spa- 
tially resolved images of the 21cm emission; besides the 
HI mass, this also enables us to determine the velocity 
field of the galaxy and, hence, to estimate its dynamical 
mass. We obtain a total HI mass of ^ 1.4 x 10^ M0, con- 
siderably smaller than the value of MJIjj determined from 
blind 21cm emission surveys. This continues the trend of 
low HI mass in all low redshift DLAs for which HI emis- 
sion observations have been attempted. We also find that 
the QSO lies behind a region of low local HI column den- 
sity in the foreground galaxy. This is interesting in view 
of suggestions that DLA samples are biased against high 
HI column density systems. Finally, the dynamical mass 
of the galaxy is found to be Mdyn 5 x 10^ M©. 

Key words: galaxies: evolution: - galaxies: formation: 
- galaxies: ISM - cosmology: observations - radio lines: 
galaxies 



1. Introduction 

Neutral gas at high redshifts is easiest to detect in absorp- 
tion against bright background sources. Not surprisingly, 
most of what we know about the content and evolution of 
neutral gas in the universe comes from the study of absorp- 
tion lines seen in the spectra of distant QSOs. The num- 
ber density (per unit redshift) of these absorption lines is 
a strong function of the column density, with low column 
density (Ahi ~ 10^"^ atoms cm^^) systems being several 



orders of magnitude more common than high column den- 
sity (Ahi ^ 10^*^ atoms cm~^) systems. Nonetheless, the 
bulk of the neutral gas at high redshift is contained in 
these rare high HI column density absorbers. It is princi- 
pally for this reason that these objects (called damped 
Lyman-a systems or DLAs) are obvious candidates for 
the precursors of today's galaxies. Further, the gas mass 
DLAs at z 3 is comparable to the stellar mass 
(e.g. Storrie-Lombardi et al. 1996 



m UIjAs at z ' 
in galaxies at z 



ptorrie-Lombardi fc Wolfe 2000] ), consistent with the idea 
that the absorbers have converted their gas into stars in 
the intervening period. Understanding the nature of DLAs 
at different redshifts is clearly important in tracing the 
evolution of galaxies; for this reason, the absorbers have 
been the subject of considerable study over the last two 
decades. 

Unfortunately, since QSOs are point sources, opti- 
cal absorption studies alone are unable to constrain the 
transverse size of the absorbers. The typical size and 
mass of DLAs have hence long been controversial issues. 
Traditionally, DLAs have been modelled as large proto- 
spirals (Wolfe et al. 1986). Some support for this model 
comes from the shapes and widths of the absorption pro- 
files produced by ions such as Sill (which are associ- 
ated with neutral HI). The large velocity widths {AV ~ 
300 km s~^) and the asymmetric shapes of these lines 
have been successfully modelled as arising from gas in a 



thick spinning disk ( [Prochaska fc Wolfe 1997| , 1998). How- 
ever, models involving infall or random motions of smaller 
gas clouds have also been found to succesfuUy repro- 



duce the observed velocity profiles (Haehnclt et al. 1998 
McDonald fc Miralda-Escude 1999]). 



At low redshifts, the galaxies in which the damped 
absorption arises can be directly studied by ground- 
based or HST observations. Contrary to expecta- 
tions, low redshift DLAs appear to be associated 
with a wide variety of galaxy types, including 
dwarf and low surface brightness (LSB) galaxies (e.g. 
Le Brun et al. 1997, Nestor et al. 2001) and not exclu- 
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sively (or even predominantly) with spiral galaxies. 
Besides this, the majority of DLAs tend to have 
low mctallicities 0.1 solar) at all redshifts, with 
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very little evolution in their metallicity with redshift 



(Pettini et al. 1999); further, they also do not show the 
expected a/Fe enrichment pattern expected for spi- 
ral galaxies ( penturion ct al. 2000| ), suggesting a differ- 
ent star formation history than that of spirals. Finally, 
21cm absorption studies ( Chengalur fc Kanekar 2d00| , 



Kanekar & Chengalur 2001) have shown that the major 
ity of DLAs have high spin temperatures (Tg ^ 1000 K), 
far higher than those observed in the Milky Way or lo- 
cal spirals (Tg ~ 200 K). All these results indicate that 
damped absorption is likely to originate in all types of 
galaxies and not merely in luminous disk systems. 

Of course, the above results do not preclude the pos- 
sibility that the galaxies responsible for the damped ab- 
sorption are indeed massive gas-rich disks, but have not 
undergone much star formation and hence have both low 
luminosities as well as low metallicities. Such low surface 
brightness systems could well have a larger fraction of the 
warm phase of neutral HI and hence, a high spin tem- 
perature ( [Chengalur fc Kanekar 2000 ). However, for sys- 
tems at low redshift, 21cm emission studies can be used 
to get direct estimates of the HI mass of the absorbers; 
one can thus directly test the above hypothesis, that 
DLAs are massive, low luminosity galaxies. Deep searches 
for HI emission from two nearby DLAs have resulted in 
non-detections (Kanekar et al. 2001, Lane 2000| ). In both 
cases, the 3 cr upper limit to the HI mass is ^ 1/3 the 
HI mass of an L* spiral. Thus, in these two cases at least, 
the absorption does not arise in an optically faint, but 
extremely gas rich galaxy. 

In this paper, we discuss Giant Metrewave Radio 
Telescope (GMRT) 21cm emission observations of a 
third low redshift DLA, the z = 0.009 absorber to- 
wards the QSO HS 1543-1-5921. The damped absorp- 
tion has been identified as arising in a low surface 
brightness galaxy SB S 1543+593 ( [Reimers fc Hagen 1998 , 
Bowen et al. 2001a). HI emission has also been de- 
tected from this galaxy using the Bonn Telescope 



( Bowen et al. 2001bD . The extremely low redshift of the 
absorber allows us to make, for the first time, spatially re- 
solved images of the 21cm emission; besides the HI mass, 
this also enables us to determine the velocity field of the 
galaxy and, hence, to estimate its dynamical mass. 

The GMRT observations are detailed in Sect. 0, while 



the results are presented in Sect. 3.1 and discussed in 



Sect. 5.2 . Throughout the paper, we use a Hubble constant 
ofHo== 75 kms^i Mpc-\ i.e. a distance to SBS 1543+593 
of 38 Mpc. 

2. Observations 

SBS 1543+593 was observed on 31 December 2000 us- 
ing the GMRT. The observing setup gave a total of 128 
spectral channels over the total bandwidth of 2 MHz, cor- 
responding to a channel separation of ~ 3.1 km s~^. The 
standard calibrators 3C48 and 3C286 were observed at the 



start and end of the observing run and used to calibrate 
the visibility amplitudes and the bandpass shape. Phase 
calibration was done using observations of the nearby con- 
tinuum sources 1438+621 and 1549+506, one of which was 
observed once every 35 minutes. 

The data were converted to FITS from the raw tele- 
scope format and analysed in classic AIPS. The brightest 
continuum source in the field of view is the steep spectrum 
source 7C 1543+5920. The flux as measured at GMRT is 
42 mjy , in excellent agreem ent with that obtained in the 
NVSS ( Condon et al. 1998 ). The source is at a distance of 
2.2 (corresponding to an impact parameter of ^ 24.1 kpc) 
from SBS 1543+593. No statistically significant 21cm ab- 
sorption is seen against this source; the 3 a limit to the 
optical depth after smoothing to a velocity resolution of 
10 km s~^ is ~ 0.1. This yields column density limits of 
~ 1.8 X lO^*' cm^^ and ^ 1.8 x 10^^ cm~^ for HI with spin 
temperatures of 100 K and 1000 K respectively, along this 
line of sight. As an aside, we note that 7C 1543+5920 has 
a spectral index of ~ —0.5 and is, at best, marginally re- 
solved on our longest baselines - the upper limit on its 
angular size is ~ 1 . Given its relatively high galactic lat- 
itude (~ 46°), it is hkely that 7C 1543+5920 is a compact 
steep spectrum (CSS) source. 

The GMRT has a hybrid configuration; 14 of its 30 
antennas are located in a central compact array with size 
^ 1 km ('^ 5 kA at 21cm) while the remaining antennas 
are distributed in the three arms of a Y configuration, 
giving a maximum baseline of ^ 25 km (^ 120 kA at 
21cm). The baselines obtained from antennas in the cen- 
tral square are similar in length to those of the VLA in its 
D configuration, while the arm antennas provide baselines 
similar in length to the B array of the VLA. A single obser- 
vation with the GMRT hence yields information on both 
large and small angular scales^]. The search for absorption 
against the compact source 7C 1543+5920 was made us- 
ing all available baselines. However, the HI emission from 
SBS 1543+593 is heavily resolved on the longest baselines. 
To study this emission, cubes were hence made at a vari- 
ety of UV ranges, viz. 0-5 kA, 0-10 kA and 0-20 kA, corre- 
sponding to spatial resolutions of 39 x 39 , 29 x 25" and 
13 X 11 . The two low resolution cubes were CLEANed, 
using the AIPS task IMAGR. At the highest resolution, 
the signal to noise ratio is too low and the emission too 
diffuse for CLEAN to work reliably. However, at the signal 
to noise ratio of the map at this resolution, convolution 
by the dirty beam does not greatly degrade the dynamic 
range (or image fidelity) . The morphology should hence be 
accurately traced, apart from an uncertainty in the scaling 
factor (which is discussed in more detail below). 

3. Results and Discussion 



^ Note, however, that the VLA in its D array has a larger 
number of short spacings than the central compact array of 
the GMRT. 
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Fig. 1. Integrated HI emission maps of SBS1543+593. The QSO position is marked by a cross in all panels. See Sect. |^ 
for a discussion of the conversion to atoms cm~^. [A] 39 x 39 resolution. The contours are 0.03, 0.045, 0.06,0.09, 
0.12, 0.18, 0.24, 0.36, 0.48, 0.72 and 0.85 Jy Beam~^ km s^^.A 3(t feature, 2 channels wide would have an integrated 
flux of 0.07 Jy Beam-i km s^^p] 29" x 25" resolution. The contours are 0.045, 0.06, 0.09, 0.12, 0.18, 0.24, 0.36 and 
0.45 Jy Beam~^ km s~^ . A 3cr feature, 2 channels wide would have an integrated flux of 0.06 Jy Beam~^ km s"'^ [C] 12 x 
11 resolution. The contours are 0.03, 0.045, 0.06, 0.09, 0.12 and 0.22 Jy Beam~^ km s~^. A 3a feature, 2 channels wide 
would have an integrated flux of 0.045 Jy Beam^^ km s~^[D] Integrated HI 21cm emission profile for SBS 1543+593. 
The profile is derived from the GMRT observations, using the 29 x 25 resolution data. The spectrum has been 
smoothed to a resolution of ^ 10 km s^^. 



3.1. Results a cross in all three maps. The centre of the optical galaxy 

is very close to the QSO position ; the QSO is located onl y 
Fig. g shows maps of the integrated HI emission (obtained 2.4" NNE of the galaxy centre ( |Reimers fc Hagen 1998D . 
using the AIPS task MOMNT) at the three different res- ^he low resolution map shows that the HI distribution 
olutions discussed above. The QSO position is marked by 
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is asymmetric with the peak emission displaced to the 
south-east of the galaxy centre. The higher resolution 
maps show clearly that the HI distribution is concen- 
trated in a ring-like structure, with HI emission actually 
being depressed at the centre of the galaxy. The HI ring 
is coincident with the faint spiral arms seen in the opti- 



cal images (Reimers & Hagen 1998, Bowen et al. 2001a). 
Many of the irregular patches of emission seen in the 
HST image ( p3owen et al. 2001a ) also appear to be asso- 
ciated with peaks of HI emission. The HI concentration 
at 15''44"20" , 59''l2'09" corresponds to the HII region 
whose spectrum is given in Reimers & Hagen (1998). The 
velocity we measure at the location of this HII region is 
2855 ± 6 km s~^. In addition to the inner HI ring, there 
are also spurs in the HI emission (see Figs. Qb] & [C]) to- 
wards the north and south; these may mark the beginning 
of faint outer spiral arms. 

Measurement of the integrated flux corresponding 
to weak extended line emission can be a non-trivial 



problem (see, for example, Jorsater & van Moorsel 1995, 
Rupen 1999 ). However, as a rule of thumb, (apart from 



serious zero spacing problems) , deeply cleaned images give 
a fairly reliable estimate of the total flux. The integrated 
fluxes that we get from the cleaned 39 x 39 and 29 x 25 
data cubes are 4.0±0.4 Jy km s~^ and 3.6±0.4 Jy km s~^. 
These are in excellent agreement with the single dish mea- 
surement of 4.0 ± 0.4 Jy km s"^ ( |Bowen et al. 2001bD . 
From the area of the clean beam in these two low resolu- 
tion images, the conversion from 1 Jy Beam~^ km s~^ to 
atoms cm~^ is 7.3 x 10^" atoms cm^^ and 1.5 x 
10^"'^ atoms cm~^ respectively. Since the highest resolu- 
tion map (Fig. |l|[C] ) is uncleaned, there is no obvious way 
to convert the units from Jy Beam~^ to Jy. If one simply 
uses the area of the best fit Gaussian to the main lobe 
of the dirty beam to convert from Jy Beam~^ to Jy, one 
gets a total flux of ~ 4.4 Jy km s^^. This must be an 
overestimate, since it is clear from Fig. |l|[C] that a con- 
siderable fraction of the smooth emission seen in Fig. ^[A] 
has been resolved out but the estimated flux in the higher 
resolution image is higher than that in the lower resolution 
ones. Given this problem in scaling for the highest reso- 
lution image, we estimate the column density at the QSO 
position only from the two lower resolution maps. The col- 
umn densities are 5.9 x 10^° atoms cm^^ and 4.9 x 10^'' 
atoms cm^^ for the 39 x 39 and the 29 x 25 res- 
olutions, respectively. The integrated emission profile of 
SBS 1543-1-593, as derived from the 29" x 25" resolution 
data (smoothed to a velocity resolution of 10 km s^^), 
is shown in Fig. |^[D]. The systemic velocity as measured 
from the profile is 2862± 10 km s^^. This agrees (within 
the errors of the two measurements) with the velocity of 
2868± 2 km s^^measured from the Bonn spectrum by 



Bowen et al. 2001a 



ric; they are straight on the approaching side and curved 
on the receding side. This type of velocity field has been 
dubbed 'kinematically lopsided' ( ^ waters et al. 1999| ). As 
noted earlier, the galaxy is also morphologically lopsided. 
The rotation curve was derived separately for the ap- 
proaching and receding sides using the AIPS task GAL. 
During these fits, the galaxy centre was kept fixed at the 
optical centre. The systemic velocity was kept fixed at the 
value of 2870 km s~^, the inclination at 50° and the posi- 
tion angle (of the receding half of the major axis, measured 
east of north) at —16° (all of which were obtained from 
an initial global fit to the velocity field). The inclination 
and the position angle are in good agreement with our 
estimates from the optical image presented in Reimers & 
Hagen (1998). 

The derived rotation curves are shown in Fig. ^[B]. 
As expected from the iso-velocity contours, the rotation 
curve is almost linear on the approaching side, while it 
tends to flatten out on the receding side. The maximum 
(inclination corrected) rotation speed in SBS 1543-f 593 is 
■~ 45 km s^^. The rotation curve can be measured out 
to a radius of ~ 60 , corresponding to a linear distance 
of 11 kpc. The implied dynamical mass is then Mdyn ~ 
5 X 10^ Mq. 

3.2. Discussion 

The flux that we measure for SBS 1543-1-593 corresponds 
to a total HI mass of 1.4 ± 0.14 x lO^M©, in excellent 
agreement with the single dish measurement of Bowen et 
al. (2001b). This is thus the third DLA whose HI mass is 
significantly less than that of an spiral. Besides these, 
the galaxy NGC 4203 (which lies in front of the QSO 
Ton 1480, and is fikely to be a DLA; [Miller et al. 1999D 



The velocity field of SBS 1543-^593 (derived using the 
AIPS task MOMNT on the 29" x 25" resolution cube) is 
shown in Fig. |^[A] . The iso-velocity contours are asymmet- 



has also been found to have a low HI mass. Thus, all low 
redshift DLAs (or candidate DLAs) for which HI emission 
observations have been attempted have masses less than 
that of an L* spiral. 

The derived HI masses are also considerably less than 
the MJjj (~ 7 X 10^ Mq) obtained from Schecter function 
fits to the HI mass function determined in blind HI surveys 
(Zwaan et al. 1997, Rosenberg & Schneider 2001). It also 
follows from these mass functions that large galaxies make 
the major contribution to the local HI mass density. How- 
ever, in order to determine the contribution of galaxies of 
different HI masses to the cross-section for DLA absorp- 
tion, one also needs to know the typical sizes of their HI 
disks. Rao & Turnshek (1998) deduced, based on determi- 
nations of the typical sizes of HI disks made by Rao (1994), 
that bright spiral galaxies make the major contribution 
to the z — cross-section for DLA absorption. However, 
more recent determinations of the HI size and mass distri- 
bution of galaxies indicate that the cross-section for DLA 
abs orption is, in fact, not domina t ed by large spiral g alax- 
ies ( Rosenberg fc Schneider 2001 , Zwaan ct al. 2001 ) , and 
that a sample of low redshift DLAs should contain a large 
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Fig. 2. [A] The velocity field of SBS 1543+593 derived from the 29 x 25 resolution cube. The velocity contours 
go from 2830 km to 2900 km and are spaced 5 km s~^ apart. Note that the contours are straight on the 
approaching (lower velocity) side and curved on the receding (higher velocity) side. The QSO position (which is only 
2.25 from the optical centre of the galaxy) is marked with a cross. [B] The rotation curve derived from the velocity 
field shown in [A]. The curve for the approaching side is marked by stars, and for the receding side by open pentagons. 
The rotation curve for the approaching side can be seen to rise approximately linearly, while that for the receding side 
fiattens out. 



variety of galaxy types. This is consistent with the optical 
and HI observations of low z DLAs. 

It would be interesting to compare the velocity profiles 
of low ionization metal lines in the HS 1543+5921 spec- 
trum with the large scale kinematics of SBS 1543+593. 
Unfortunately, high velocity resolution optical/UV ob- 
servations are as yet not available for this system. 
We note, however, that the systemic velocity that we 
obtain for SBS 1543+593 is in excellent agreement 
with the value of 2868 ± 2 km s^^ obtained from 
the single dish HI profile ( Bowen et al. 200rh| ). Fur- 
ther, the velocity obtained in our observations (2855 ± 
6 km s^^) for the gas coincident with the HII region 
(for which an optical emission spectrum exists) is in rea- 
sonable agreement with that obtained in the optical (~ 
2700 km s~^; Reimers & Hagen 1998), given the poor res- 
olution (18 angstroms, i.e. ~ 1000 km s~^) of the op- 
tical spectrum. Similarly, the velocity we obtain at the 
QSO location is ~ 2870 ± 10 km s^^, which agrees within 
the error bars with the velocity of 2700 km s~^ ob- 



kincmatics of the absorbing galaxies (see, for example. 



tained from the Ly— a line (Bowen et al. 2001a), given 
the large uncertainty (~ 200 — 300 km s^^) in the lat- 
ter measurement (Bowen et al. 2001b). A high resolu- 



tion low ionization metal line absorption spectrum of the 
QSO (and of Ton 1480, which lies behind NGC 4203) 
would provide interesting spot checks of the reliability 
in using such absorption lines to probe large scale gas 



Prochaska fc Wolfe 1997, Prochaska fc Wolfe 199S) 



It is also interesting to compare the HI column density 
as obtained from the Ly-a spectrum with that obtained 
from the 21cm emission. The comparison is, however, com- 
plicated by the comparatively large size of the GMRT syn- 
thesized beam. At a distance of 38 Mpc, 29 corresponds 
to a linear size of ~ 5.3 kpc. Given this great disparity 
in the transverse sizes covered by the HI and UV mea- 
surements, the column density of 5 x 10^° atoms cm~^ 
obtained from the 21cm map is in reasonable agreement 
with the estimate of 2.2 x 10^" atoms cm~^ obtained from 
the Lyman-a profile. 

Comparisons of DLA samples with samples of galax- 
ies selected by blind HI emission are complicated by 
issues of optical extinction. It has long been sug- 
gested that DLA samples could be seriously biased 
against high HI column density systems, because such 
absorbers would contain enough dust to substantially 



dim the background QSO (e.g. Hcislcr & Ostriker 1988 
Pei et al. 1999). Deep optical observa- 



FaU & Pei 1993 



tions of a radio selected sample of QSOs however suggest 



that such biases may be modest (Ellison et al. 2001). It is 



curious (but, of course, not statistically significant, partic- 
ularly given the fact that optical depth effects may be less 
important in LSB galaxies than in normal spirals) that 
the QSO HS 1543+5921 lies behind a region of low lo- 
cal HI column density in SBS 1543+593. Interestingly, an 
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inspection of the HI maps of the SO galaxy NGC 4203 pre- 
sented by van Driel et al. (1988) shows that Ton 1480 also 
lies behind a region of low column density. Unfortunately, 
testing the above bias would require 21cm mapping of a 
statistically significant number of DLAs which is, alas, a 
task for the next generation radio telescope. 
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